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^ ' Abstract. We study the Cauchy problem for a scalar semilinear degenerate 

parabolic partial differential equation with stochastic forcing. In particular, 
we are concerned with the well-posedness in any space dimension. We adapt 
the notion of kinetic solution which is well suited for degenerate parabolic 
problems and supplies a good technical framework to prove the comparison 
principle. The proof of existence is based on the vanishing viscosity method: 
^■f^ , the solution is obtained by a compactness argument as the limit of solutions 

■ of nondegenerate approximations. 
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j_j 1. Introduction 



In this paper, we study the Cauchy problem for a scalar semilinear degenerate 
parabolic partial differential equation with stochastic forcing 

du + div{B{u))dt = diy{A{x)\7u)dt + ${u)dW, a; G T^, t G (0, T), 

■ti(O) = Uq, 

where is a cylindrical Wiener process. Equations of this type are widely used 
in fluid mechanics since they model the phenomenon of convection-diffusion of 
ideal fluid in porous media. Namely, the important applications including for 
instance two or three-phase flows can be found in petroleum engineering or in 
hydrogeology. For a thorough exposition of this area given from a practical point 
of view we refer the reader to lUj and to the references cited therein. 

The aim of the present paper is to establish the well-posedness theory for 
solutions of the Cauchy problem (jl.ip in any space dimension. Towards this 
end, we adapt the notion of kinetic formulation and kinetic solution which has 



This research was supported in part by the GACR Grant no. P201/10/0752. 



2 



Martina Hofmanova 



already been studied in the case of hyperbolic scalar conservation laws in both 
deterministic (see e.g. [TH], [TH], [20], [55], or [33] for a general presentation) 
and stochastic setting (see [S]); and also in the case of deterministic degenerate 
parabolic equations of second-order (see [3]). To the best of our knowledge, in 
the degenerate case, stochastic equations of type ()l.ip have not been studied 
yet, neither by means of kinetic formulation nor by any other approach. 

The concept of kinetic solution was first introduced by Lions, Perthame, 
Tadmor in [2D] for deterministic first-order scalar conservation laws and applies 
to more general situations than the one of entropy solution as considered for 
example in [2], [S], [T7]. Moreover, it appears to be better suited particularly for 
degenerate parabolic problems since it allows us to keep the precise structure 
of the parabolic dissipative measure, whereas in the case of entropy solution 
part of this information is lost and has to be recovered at some stage. This 
technique also supplies a good technical framework to prove the L^-comparison 
principle which allows to prove uniqueness. Nevertheless, kinetic formulation 
can be derived only for smooth solutions hence the classical result giving 
L^'-valued solutions for the nondegenerate case has to be improved (see [15]). 

In the case of hyperbolic scalar conservation laws, Debussche and Vovelle 
[5] defined a notion of generalized kinetic solution and obtained a comparison 
result showing that any generalized kinetic solution is actually a kinetic solution. 
Accordingly, the proof of existence simplified since only weak convergence of 
approximate viscous solutions was necessary. 

The situation is quite different in the case of parabolic scalar conservation 
laws. Indeed, due to the parabolic term, the approach of [S] is not applicable: the 
comparison principle can be proved only for kinetic solutions (not generalized 
ones) and therefore strong convergence of approximate solutions is needed in 
order to prove the existence. Moreover, the proof of the comparison principle 
itself is much more delicate then in the hyperbolic case. 

The exposition is organised as follows. In Section 2 we review the basic set- 
ting and define the notion of kinetic solution. Section 3 is devoted to the proof 
of uniqueness. We first establish a technical Proposition 13.21 which then turns 
out to be the keystone in the proof of comparison principle in Theorem l3.3l We 
next turn to the proof of existence in Sections 4 and 5. First of all in Section 
4, we make an additional hypotheses upon the initial condition and employ the 
vanishing viscosity method. In particular, we study certain nondegenerate prob- 
lems and establish suitable uniform estimates for the corresponding sequence 
of approximate solutions. The compactness argument then yields the existence 
of a martingale kinetic solution which together with the pathwise uniqueness 
gives the desired kinetic solution (defined on the original stochastic basis). In 
the final section, the existence of a kinetic solution is shown for general initial 
data. 

We note that an important step in the proof of existence, identification of 
the limit of an approximating sequence of solutions, is based on a new general 
method of constructing martingale solutions of SPDEs (see Propositions 14.141 
14. isl and the sequel), that does not rely on any kind of martingale representation 
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theoreni and therefore holds independent interest especially in situations where 
these representation theorems are no longer available. First applications were 
already done in [T], [21] and, in the finite-dimensional case, also in jH] , 

2. Notation and main result 

We now give the precise assumptions on each of the terms appearing in the above 
equation (II. We work on a finite-time interval [0,r], T > 0, and consider 
periodic boundary conditions: x G where is the A^-dimensional torus. 
The flux function 

B = {Bi,...,Bn) :M — 

is supposed to be of class with a polynomial growth of its derivative, which 
is denoted hy b = (6i, . . . , b^). The diffusion matrix 

is of class C°°, symmetric and positive semidefinite. Its square- root matrix, 
which is also symmetric and positive semidefinite, is denoted by a. 

Regarding the stochastic term, let (fi, {^t)t>o,^) be a stochastic basis 
with a complete, right-continuous filtration. The initial datum may be random 
in general, namely, we assume uq S Lp{^; LP{T^)) for all p £ [l,c)o). The 
process is a cylindrical Wiener process: W{t) — J2k>i f^k{t)ek with {I3k)k>i 
being mutually independent real-valued standard Wiener processes relative to 
{•^t)t>o and {ek)k>i a complete orthonormal system in a separable Hilbert space 
il. In this setting, we can assume, without loss of generality, that the cr-algebra 
^ is countably generated and {^t)t>o is the filtration generated by the Wiener 
process and the initial condition. For each z e L^(T^) we consider a mapping 
<P{z) : il ^ L^{T^) defined by <P{z)ek = gk{-,z{-)). In particular, we suppose 
that gk G C(T^ x R) and the following conditions 

G2(x,e) = 5^|ff.(x,C)|'<i(l + |e|^), (2.1) 

k>l 

E - 9kiy,0\^ < H\x - + \^- CIMI^ - CD), (2.2) 

k>l 

are fulfilled for every a;, y G T^, ^, ^ G K, where h is a continuous nondecreasing 
function on ]R_|_ satisfying, for some a > 0, 

h{6) < CS"', 6<1. (2.3) 

The conditions imposed on (p, particularly assumption (|2.ip . imply that 

where L2{ii; L'^{T^)) denotes the collection of Hilbert-Schmidt operators from 
il to X^(T^). Thus, given a predictable process 

ueL\n;L^O,T;L\T''))), 
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the stochastic integral t M- J* <P{u)dW is a well defined process taking values in 
^^(T^) (see [5] for detailed construction). 

Finally, define the auxiliary space Hq D H via 



2 

*;>! fe>l 



endowed with the norm 

2 

A;>1 fc>l 

Note that the embedding it ilo is Hilbcrt-Schmidt. Moreover, trajectories of 
W are P-a.s. in C{[0,T];%) (see [5j. 

As the next step, we introduce the kinetic formulation of as well as 
the basic definitions concerning the notion of kinetic solution. 

Definition 2.1 (Kinetic measure). A mapping m from 51 to the set of nonnegative 
finite measures over x [0, T] x M is said to be a kinetic measure if 

(i) in is measurable in the following sense: for each ijj G Cf,(T^ x [0,r] x R) 
the mapping (m, "0) : J7 — >■ M is measurable, 

(ii) m vanishes for large ^: if Bf^ = € R; |^| > R} then 

lim Eto(T^ X [0,T] X B%) = 0, (2.4) 

R^oo 

(hi) for aU V e Cb(T^ x M), the process 

Jt" x[0,i]xR 

is predictable. 
Definition 2.2. The kinetic formulation of (jl.ip is 

dtf+biO-^f-Y. =^^=i^i^)^+9Am~-G^Su=A, (2.5) 

where m is a kinetic measure which consists of two components m — ni + n2. 
The measure ni is known and relates to the diffusion term in (II. ip 

ni{x,t,^) = (Vu)*A(a;)(Vu)5„=5 

whereas n2 is an unknown nonnegative measure. 

Remark 2.3. The measure ni is called parabolic dissipative measure and gives 
us better regularity of solutions in the nondegeneracy zones of the diffusion 
matrix A (cf. Definition 12.41 (iii)). The measure n2 takes account of possible 
singularities of solution and vanishes in the nondegenerate case. 

We now derive the kinetic formulation in case of a sufficiently smooth u 
satisfying It is a consequence of the Ito formula. 



Degenerate Parabolic SPDEs 



5 



Let us set (ln>{,6'') = l„>46l'(OdC = 0{u) for 9 e Cf'(M) and apply 
the Ito formula: 

d(l„>5 , 6*') 9'{u) [ - div {B{u))dt + div {A{x)\7u)dt + ${u)dW{t)] 
+ ^9"{u)G\u)dt. 
Afterwards, we proceed term by term 
9'{u) div {B{u)) = 9\u)b{u)- Vu = div ^ 6(^)6'' (^d^^ = div ((61„>4, 6*')) 

N N 

9' {u) div {A{x)Vu) = d^^[A,j{x)9'{u)d^^u] - ^ 9"{u)^^M^A^)9x,u 

^'H.Z'Cu) = {Su=^$iu),9') 

9"{u)G^u) = (G2<5„=^,0") = -{d^{GHu=d,9') 

Taking 6l(^) = for any test function ip e C^iJ^ x [0,r] x R) we 

obtain that / = lu>^ is a distributional solution to the kinetic formulation 
p.Sp with 712 = 0. Therefore any smooth solution of (jl.ip is a kinetic solution 
in the sense of the following definition. 

Definition 2.4 (Kinetic solution). A measurable function m : O x x [0, T] ^ M 

is said to be a kinetic solution to (11.11) with initial datum ug provided 

(i) t £ [0,T]) is a predictable L^(T^)-valued process, p G [l,c)o), 

(ii) for all p G [1, oo) there exists Cp > such that for t G [0, T] 

Mt)\\LPinxT^)<Cp, (2.6) 

(iii) {\7u)*A{x){\7u) G Li(17 X X [0,r]), 

(iv) there exists a kinetic measure m such that m = ni + n2 as in Definition 
land / = lu>4 is a weak solution of the kinetic formulation ()2.5p . i.e. 

-a.s. for ah (p G G^{T^ x [0,r] x R), 

{f{t),dMt))dt + (/o, ^(0)) + r (fit), b{0- V^{t))dt 

Jq 

+ / (/W, E 5.,(A,,(x)a.,^(t)))di 

"'O i,j=l 

E / f ^gk{x,u{x,t))ip{x,t,u{x,t))dxdl3kit) 

G^(x, w(a;, i))9^(/3(x, i, u(a;, t))dx dt + m{d^ip). 



k>l 



/o JT" 



'0 "'T" 

(2.7) 
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We proceed by two related definitions, which will be useful especially in 
the proof of uniqueness. 

Definition 2.5 (Young measm-e). Let {X, A) be a finite measure space. A map- 
ping 1/ from X to the set of probability measures on K is said to be a Young 
measure if, for all ip £ Cf,(K), the map z i— >■ :^z(V') from X into R is measurable. 
We say that a Young measure ly vanishes at infinity if, for all p > 1 , 



X 



1^1^1^,(0 dA(z) < 



oo. 



Definition 2.6 (Kinetic function). Let {X, A) be a finite measure space. A mea- 
surable function / : X x M — > [0, 1] is said to be a kinetic function if there exists 
a Young measure on X vanishing at infinity such that, for A-a.e. z d X, for 
all C G R, 

Remark 2.7. Note, that if / is a kinetic function then d^f — --v. Similarly, let 
u be a kinetic solution of (jl.ip and consider / — lu>j. We have d^f — —6u=^, 
where — Su=^ is a Young measure on x x [0, T]. Therefore, the expression 
((T7| can be rewritten in the following form: for ah ip G C;?°(T^ x [0,T] x R), 
P-a.s., 

^{fit),dMt))dt + {fo,m) + r {fit)MO- v^(0)di 

Jo 

N 



+ I {.fit), E d.AA^ii^)9.Mt)))dt 
[ [ [ 9k{x,0^{x,t,Odiy.,^t{Odxdl3k{t) 



T 

{x, ^)d^(p{x, t, Odi^x.t (C) dx dt + m{d^ip). 

Jt" 



(2.8) 

For a general kinetic function / with corresponding Young measure v, the above 
formulation leads to the notion of generalized kinetic solution as used in j5] . Al- 
though this concept is not established here, the notation will be used throughout 
the paper, i.e. we will often write Vx.tiO instead of 5ti(a;,t)=^- 

Lemma 2.8. Let (X, A) be a finite measure space such that L^{X) is separable^ 
Let {/„; n G N} he a sequence of kinetic functions on X x R, i.e. fn{z,C) = 
t'"(^, 00) where are Young measures on X. Suppose that, for some p > 1, 



sup / / \^\Pdv^{Od\{z)<^. 

neNJx JR 



^According to [4] Proposition 3.4.5], it is sufficient to assume that the corresponding cr-algebra 
is countably generated. 
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Then there exists a kinetic function f on X xM. and a subsequence still denoted 
by {fn', n Cz N} such that 

fn ^ f, m L°°{X X R)-weak*. 

Proof. The proof can be found in 6, Corollary 6]. □ 

To conclude this section we state the main result of the paper. 

Theorem 2.9. Let uq e LP(r2; LP(T^)), for all p e [l,oo). Under the above as- 
sumptions, there exists a unique kinetic solution to the problem (II. ip . Moreover, 
ifui, Ui are kinetic solutions to (jl.ip with initial dataui^ andu2fl, respectively, 
then for all t G [0,r] 

El|ui(i) - M2(i)llLi(T") < El|ui,o - U2,o1Il1(T")- 

3. Uniqueness 

We begin with the question of uniqueness. Due to the following proposition, we 
obtain an auxiliary property of kinetic solutions, which will be useful later on 
in the proof of the comparison principle in Theorem 13.31 Since the proof is very 
similar to [SJ Proposition 8], it will be left to the reader. 

Proposition 3.1 (Left and right weak limits). Let u he a kinetic solution to (jl.ip . 
Then f = 1m>^ admits almost surely left and right limits at all points t* G [0, T] 
in the sense of distributions over X R, i.e. for all t* G [0, T] there exist some 
kinetic functions f*'^ on fl x x R such that 

(/(r±e),^)^(r'±,V), e;0, VV-GCriT^xR) P-a.s.. 

Moreover, f*'^ = f*'^ = f{t*) almost surely for all t* G [0,T] except for some 
countable set, i.e. 

(r'+,V') = (r^~,^) = (/(**), v-) vv^Gcr(T^xR) p-a.s. 

As the next step towards the proof of the comparison principle, we need 
a technical proposition relating two kinetic solutions of (II. ip . We will also use 
the following notation: if / : X x R — > [0, 1] is a kinetic function, we denote by 
/ the conjugate function / = 1 — /. We define the function by f^{t*) = 
t* G [0,T]. 

Proposition 3.2. Let ui, U2 be two kinetic solutions to (jl.ip and denote fi ~ 
ltii>Ci I2 — lii2>4- Then for t G [0,r] and any nonnegative functions g G 
C°°(T^), V G C^{M.) we have 

e[ f g{x~y)iji^~C)fti^,t>0f2{y>tX)d^dCdxdy 

<e/ / 0/1,0 (a^, 0/2.0 (2/,C)dCdCdxdy + I + J + K, 

(3.1) 
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where 



l^^f f I /i/2(fe(0-KC))-V,a(x,e,y,C)dedCdxd2/ds, 



J=E [ f I hhiZ^v. (^^J ^y)^v. «) dC da; dy As 

Jo i(T«)2 JR2 . 

+ E / / / /1/2 V 9,,.(A,,(x)5,.a)dedCdxdyds 

f f I a(x,e,y,C)d^^^'+(0da;dn2,i(y,s,C) 
-E / a(a;,^,y,C)di^^'7(C)dydni.i(x, s,^), 

■/O J(T")2 JR2 



K^Ie/Y / a(x,e,y,C)Ek'=(^'0-5fc(y,C)|'d^^.«(0d^^L(C)d^dyds, 



anc? the function a is defined as a{x, ^, y, (^) — g{x — y)ijj{S, — C)- 



Proof. Let us denote by ((•, •}} the scalar product in ^^(Tf x x x R^). 
In order to prove the statement in the case of f^, f2, we employ similar calcu- 
lations as in Proposition 9] to obtain 
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+ e/ / / /i/2(5(0"6(C))-V,ad^dCdxd2/ds 
+ E / / I hhy\ dy^ {A,, {y)dy^ a) d^ dC dx dy ds 

Jo J(T«)2 JR2 . 

+ e/ / [ fj2y d^^{A,,{x)d,,a)d^d(:dxdyds 

Jo J(T«)2 JR2 

+ ^Ie/ / / f2d^aGldiylJ0dCdydxds (3.2) 

^ Jo J(T")2 JR2 

I I hdcaGldvl,{(:)d£,dydxds 

^ Jo ■/(T")2 JR2 
Jo J(T")2 JR2 

-E / /2~94admi(a;, s,^)dCdy 

Jo J(T")2 JR2 

+ E f f f f+d^adm2{y,s,C)d^dx. 

Jo J(T")2 JR2 



In particular, since a > 0, the last term in p.2p satisfies 

IE / / _ f f^dcadm2iy,s,C)d^dx 
ft 



J(T«)2 



-E / / / adj.i;+(0dxdn2,i(y,s,C) 

Jo J(T"~|2 Jr2 



/O J(T")2 JR2 

f f f adiyl:t{0dxdn2,2{y,s,0 

Jo J(T")2 JK2 

<-^l I I advl^+{£,)dxdn2Av.sX) 
Jo J(T^y Jr2 



and by symmetry 

~E / f^d^admi{x,s,()d(dy 

Jo JfT")2 Jr2 



J(T")2 JR2 

t 

,2, 



<-E/ / / ad<7(C)dydni.i(a;,s,e). 

Jo J(T")2 . 



Thus, the desired estimate p.ip follows. 

In the case of f^, f2 we take t„ 1 write (13. ip for f^{tn), f2{tn) and let 
n — > oo. □ 
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Theorem 3.3 (Comparison principle). Let u be a kinetic solution to (jl.ip . Then 
u has left and right limits at any point in the sense of L^(T^), p e [1, oo), and, 
for aUtG[0,T],f^{x,t,0 - lM±(x,t)>^ fl-S-j for a.e. (a;,^). Moreover, if Ui,U2 
are kinetic solutions to (jl.ll) with initial data ui^ and U2fi, respectively, then 
for all te[0,T] 

E\\uf{t) ~ uf{t)\\LHT«) < E||ui,o - M2,o||li(t")- (3.3) 

Proof. Denote fi = l„i>^, /2 = lti2>?- Let (gr), i^s) be approximations to 
tlie identity on and R, respectively, i.e. let G C°°(T^), V e C^{R) be 
nonnegative symmetric functions satisfying Jrj.jv £* = Ij /r^ = 1 supp C 
B{0, 1/2), suppV' C (-1, 1). We define 

Then 



E / / ft{x,t,0f2i^.t,0<i^dx 

= e/ / gr{x-yyjs{i-0fH^,t^0f2{y,tX)d^d(:dxdy + Vt{r,S), 

where limT-_5_yo ?7t (t, (5) — 0. With regard to Proposition 13.21 we need to find 
suitable bounds for terms I, J, K. 

Since b has at most polynomial growth, there exist C > 0, p > 1 such that 

H)-m\<mcM-c\, mc)<c{i + \e-' + \cr')- 

Hence 

\l\<Et f I fif2r{tO\^-Cm^-Od^dC\V,Qr{x-y)\Ax<lyds. 

Jo J(T")2 JR2 

As the next step we apply integration by parts with respect to ^. Focusing 
only on the relevant integrals we get 



/ f2{c)mx)\^-c\M^-odcdc 

= / fliO l^mx'M-C\MC-C)dCd^ 

r mO\i-C'\M^-C')dC'd^d,yl,iO (3.4) 

J — oo 

/oo 
m,m - c\M^ - c)dc d^di^l^ic) 



where 



/ m',a\e-c\Me-c)dCde. 
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\\\<^f I I r{i,Qdvl,{i)dvl,{Q\V.Qr{x-y)\dxdyds. 

Jo J(TW)2 Jr2 

The function T can be estimated using the substitution ^" = ^' — (' 

mo = r [ r{C + C,C)\e\MC) d^dC 

J( J\e'\<s,e'<i-c 

< C6 / max r(£," + C) dC 

< cs (1 + ler^ + icr-') dC 
<c5(i + icr + ic'r) 

hence, since v^, vanish at infinity, 

|I| < CM / \^xQr{x)\ dx < CtSr-'^. 
We recall that /i = l„i(x,t)>i, /2 = l«2(!/,t)>C and 

The first term in J can be rewritten in the following manner using inte- 
gration by parts (and considering only relevant integrals) 

/ /2 / fida;^{Aij{x)dxiQT{x - y))dxdy 

= -/ h{y,s,C) I dxJi{x,s,C)-^ij{x)dxiQT{x-y)dxdy 

JfN JjN 

= / f2iy,s,C)dxjfi{x,s,0Aj{x)dy.gr{x - y)dxdy 
= - / dxJi{x,s,0Aij{x)dyj2{y,s,C)QT{x - y)dxdy, 

J(T«)2 

similarly 

/ /i / f2dy^{Aij{y)dy,gr{x-y))dydx 

= - dxJi{x,s,S,)Aij{y)dyj2{y,s,C)Qr{x - y)dxdy. 

Using rcgularization by convolutions, it is possible to show that the following 
holds in the sense of distributions over x M. 

dxifi = dxiUi5u-,{x,s)=i, dyja = -9yiU2<5u2(j/,s)=f • 
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Hence 



Jo J(T")2 


- ym[ui{x 


,s) - 


- U2{y, s)) dxdyds 




y)il;s[ui(x, 


s)- 


U2{y, s)) da; dy ds 


-e/ / {'^yU2)*A{y){VyU2)grix- 


y)%l)s{ui{x, 


s)- 


U2{y, s)) dxdyds 


Jo J(T")2 


- y)ips{ui{x, 




' U2{y, s)) dxdyds. 


Let us define 









^ — oo 

Then we have J = Ji + J2 + J3 with 

Ji=-E/ / (Va;Ui)V(x)CT(x)(V£ir)(a:-y)65(ui(x, s) -U2(y, s))dxdyds, 

Jo J(T")2 

J2=-E/ / (Vj,M2)*cr(y)cr(y)(Vgr)(a;-y)05(ui(x,s) -U2(y, s))dxdyds, 

Jo J(T")2 

J3 = -e/ / [|a(x)V:,ui|2 + |a(y)V;,U2|^]£.^(x-y) 

Jo J(T")2 

X V'5('"i(a;, s) - W2(y, s)) dxdyds. 

Let 

H = e/ / {\/xUi)*(T{x)a{y){\/yU2)gT{x-y)ips{ui{x,s)-U2{y.,s)) dxdyds. 
Jo J(T^Y 

In order to prove that J is nonpositive for t small enough, it is sufficient to show 
Ji = H + 0(1), J2 = H + 0(1), where o(l) ^> as t — > uniformly in 5. Indeed, 
we then obtain 

J = -E/ / \a(x)Va:Ui - a{y)VyU2?'gT{x - y) 
Jo J(T^Y 

X -05(^1 (x, s) - U2(y, s)) dxdyds + o(l) < o(l). 

We will prove the claim for Ji only since the proof for J2 is analogous. Define 

g(x,y,s) = (Va;-ui)*f7(x)85(ui(x, s) -U2(y,s)). 

Here, we make use of the assumption (iii) in Definition l2.4l Observe, that it gives 
us some regularity of the solution in the nondegeneracy zones of the diffusion 
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matrix A and hence g G 2.2(1] x x x [0, T]). It holds 



Ji = -E 



J(T")2 



E 



H = E 



-/(T^j^ 
t 



/(x, y, s) {a{x) - a{y)^ (V£ir)(a; - y) dx dy 



g{x,y,s)a{y){V Qr){x - y) dxdyds, 

-'(T")2 

(z, y, s) divj^ {a(y)Qr{x - y)^ dxdyds 



ds 



J(1N)^ 

t 

E 



l{x, y, s) div (a{y))gr{x - y) dx dy ds 
g{x,y,s)a{y){\7gr){x - y) dxdyds, 

J{Y«Y 

where divergence is applied row-wise to a matrix-valued function. Therefore, it 
is enough to show that the first terms in Ji and H have the same limit value if 
T — !> 0. For H, we obtain easily 

g{x, y, s) div(cr(?/)) gr{x - y) dx dy ds 

J(T")2 

— ^ E / / g{y,y, s) div (a (y)) dyds 
Jo Jt" 

so it remains to verify 

l{x,y,s)(^a{x) - a{y)^{Vgr){x - y) dxdyds 

E / / g{y,y,s)dvf [(j{y)) dyds. 



J(T")2 



"'T« 



We will use similar arguments as in the commutation lemma of DiPerna and 
Lions (see [71 Lemma II.l]). Let us denote by g^ the element of g and by cr* 
the i*'* row of a. Since T|Vpr|(') £ Cg2T{-) with a constant independent of r, 
we obtain the following estimate 



E 





/ 9'ix,y,s)( 


Jo Jt" 





dy ds 



< C ess sup 


a'{x')-a\y') 


^ / \g'ix,y,s)\g2r{x - y)dxdyds 




T 


Jo J(T«)2 



\x' —y' \ <T 

Note that according to [TU], [21], the square-root matrix of A is Lipschitz con- 
tinuous and therefore the essential supremum can be estimated by a constant 
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independent of r. Next 



E 



/ / W{x,y,s)\Q2r{x-y)dxdyds 

Jo J(T«)2 

< (e / / \9'{x,y,s)\'^g2T{x -y)dxdyds] 

\ Jo J(T")2 / 

X ( / g2T{x -y)dxdy] 

/ / |(Va;Wi)*o-(a;)|^ / Q2t{x -y)dydxds 
Jo Jt" Jt" 



< ( E 



So we get an estimate which is independent of t and S. It is sufficient to con- 
sider the case when 5* and ct' are smooth. The general case follows by density 
argument from the above bound. It holds 

-E r / g\x,y)(a\x) - a\y)){WQr){x - y)dxdyds 

Jo J(T")2 ^ ^ 

= WTT^ f I I 9\x,y)Da'{y + r{x-y)){x-y) 

■ (V^) -j dr da;dy ds 



= — E / / / g\y + TZ,y)Da^{y + rTz)z ■ {V q){z) dr dz dy ds 

Jo >/(TN)2 ^0 

^ ~E f f g\y,y)T>a\y)z-{WQ){z)dzdyds. 

Jo 7(T")2 

Integration by parts now yields 

/ Zid^g{z)dz = -dij, i,jG{l,...,N}, 

hence 

-E /* / g\y, y)Da' {y)z ■ {Vg) {z)dzdyds = E f [ g\y,y) div (a* {y)) dyds 
Jo J(tjv)2 Jo Jjn 

and we deduce finally that J is nonpositive. 
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The last term K is, due to (|2.2p . bounded as follows 
K<^e/ / Qr{x - y)\x - [ ^^(^ - C) di^i,.(0 dz^L(0 dxdy ds 



J(T«)2 



<7^ \x-y\^gr{x-y)dxdy+ / 



Qr{x - J/) Ax Ay 

/(T")2 ^ J(T")2 

^Y^ " + — 2 — ' 

where = sup^^g Finally, we deduce for all t 6 [0,T] 

fHx,t,0f2ix,t.0d^dx<E f f /i,o(a;, 0/2,0(2:, Od^dx 

T" JK JT" JR 

+ Ct((5T-l + 5- V2 + + TJt{T, S) + 77o(r, (5). 

Taking 5 — t^^^ and letting t ^ yields 

E / / /±(t)/±(t)ded:E < E / / /i,o/2,odeda;. 

Let us consider now ]\ — ji — /• Since /o = luo>c have the identity /o/o = 
and therefore — /^) = a.e. (a;,a;,0 and for all t. The fact that is a 
kinetic function hence implies that there exist it^*^ : x T-*^ x [0, T] R such 
that = l„±>j for almost every (a;,x,0 and all t. Furthermore, it follows 
from Proposition 13. II that u+ = = u except for a countable set of t. Since 



/ lnf>cl.±>«de=("?-4) + 
J B 

we have the comparison property 

- '"^(t))^||^i^^jvj < E||(mi,o - U2,o)^||ii(TiV)- 

It remains to show that are also left and right limits of u in the sense of 
L^(T^). But this is a consequence of the following: for p > 1 and s,t G [0, T], s > 
t 

\u{x,s) -u+{x,t)fAx ^ I l {lu(x,s}>i-K+{x.t)>e)-^'iidx 



{fix,s,^)-f+ix,t,0)^A^Ax 
and the claim follows from the weak*-convergence of f{s) f^{t) as s J, t. □ 



As a consequence of Theorem 13.31 namely from the comparison property 
, we obtain the uniqueness part of Theorem 12.91 The proof is similar to [HI 
Corollary 121. 
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Corollary 3.4 (Continuity in time). Let u : Q x x [0, T] ^ R be a kinetic 
solution to (jl.ip . Then, for all p £ [l,oo), u has almost surely continuous tra- 
jectories in LP{T^). 

4. Existence - smooth initial data 

In this section we prove the existence part of Theorem 12.91 under an additional 
assumption upon the initial condition: uq £ LP{il; C°°(T^)), for all p £ [1, oo). 
We employ the vanishing viscosity method, i.e. we approximate the equation 
(|1.1|) by certain nondegenerate problems, while using also some appropriately 
chosen approximations of 'P and B, respectively. These equations have 

smooth solutions and consequent passage to the limit gives the existence of a 
kinetic solution to the original equation. Nevertheless, the limit argument is 
quite technical and has to be done in several steps. It is based on the com- 
pactness method: the uniform energy estimates yield tightness of a sequence 
of approximate solutions and thus, on another probability space, this sequence 
converges almost surely due to the Skorokhod representation theorem. The limit 
is then shown to be a martingale kinetic solution to (|1.1[) . Combining this fact 
and the pathwise uniqueness with the the Gyongy-Krylov characterization of 
convergence in probability, we finally obtain the desired kinetic solution. 

4.1. Nondegenerate case 

Consider approximations to the identity {(pe), (ipe) 

on X M and R, respec- 
tively, and a truncation (xe) on M, i.e. we define XeiO = xi^O^ where x is a 
smooth function with bounded support satisfying < x ^ 1 and 



where x £ T^, ^ £ M. Consequently, we set B^ = [Bf, . . . , Bf^) and define the 
operator by (l>^{z)ek = g^{-,z{-)), z £ L^(T^). Clearly, the approximations 
i?^, g\ are of class C°° with a compact support therefore Lipschitz continuous. 
Moreover, the functions g| satisfy (|2.ip . (|2.2p uniformly in e and the following 
Lipschitz condition holds true 



Vx£T^ V^, CeM ^|g|(x,0-5l(a:,C)r <^e|e-CI'- (4-1) 



From (P?T|) we conclude that ^^{z) is Hilbert- Schmidt for all z £ L'^iJ^). Also 
the polynomial growth of B remains valid for B'^ and holds uniformly in e. 




The regularizations of B are then defined in the following way 




fe>i 
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Suitable approximation of the diffusion matrix A is obtained as its pertur- 
bation by el, where I denotes the identity matrix. We denote A'^ ^ A + el. 

Consider an approximation of problem (jl.ip by a nondegenerate equation 

du' +dW{B''{u''))dt^div{A^{x)Vu')dt + <P^ {u^ ) dW, 

u%0)=uo. ^'^•^^ 

Theorem 4.1. Let uq e LP{n; C°°(T^)), for all p e [1, oo). For any e > 0, there 
exists a predictable C°°{T^) -valued process solving 



Proof. For any fixed e > 0, the assumptions of [HI Theorem 2.1] are satisfied 
and therefore the claim follows. □ 

Thus, using the same computation as in the Section [H one can verify that 
the solution satisfies the kinetic formulation of (j4.2|) : let f^ = l„e>^ 

N 

where m'^ = nf + nl and both these measures are explicitly known and corre- 
spond to the diffusion matrix A + el: 

Note, that by taking limit in e we lose this precise structure of n2 . As can be seen 
from the derivation of kinetic formulation, u'^ solve ()2.5p even in a stronger sense 
than the one from Definition O Indeed, for any Lp G C'^{T^ x M), t £ [0, T], 
the following holds true P-a.s. 

(r v) - (/o, ^) - f (r (s), v^) ds 

Jo 

-/ E a,^(A,,(a;)a,,^))ds-e / {f^.s),A^)ds 

Jo \ / ^0 

{Su^=i $'{u') dW, 'P) + lf {Su'^i G\ d^v) ds - (m^ 95^)([0, t)). 

(4.3) 

4.2. Energy estimates 

In this subsection we shall establish the so-called energy estimate that makes 
it possible to find uniform bounds for approximate solutions and that will later 
on yield a solution by invoking a compactness argument. 

Lemma 4.2. For all e G (0,1), for all t e [0,T] and for all p e [2,oo), the 
solution u'^ satisfies the inequality 

Eh'(i)IIL(T«) < C(l + E||uo|l^,(T«))- (4-4) 
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Proof. We apply the Ito formula using f{v) ~ Ikll^p(Tiv)- If Q is the conjugate 
exponent to p then f'{v) — p\v\P^'^v G _L'(T^) and 

f"{v)^p{p-l)\v\P~^ld e^(LP(T^),L«(T^)). 

Therefore 

ll"^WIIL(T«) =ll"o|li.(T«) -P r / \u^P-'u^div{B^u^))dxds 
\ ' ^ ' Jo Jt« 

+ P f f \u'\P^^u'dw{A{x)Vu')dxds 
Jo Jt" 

+ ep [ [ \u^\P~'^u''Au^ dxds (4.5) 

We conclude that the second term on the right hand side vanishes and using 
the integration by parts, the third one as well as the fourth one is nonpositive. 
The last term is estimated as follows 



]-p{p-l) f f \u'\P-^Gl{x,u')dxds<C f f \u'\P-^l + \u'\^)dxds 
2 Jo Jt™ Jo Jt" 

<c(l + ^V^(s)IIL(T«)ds). 

Finally, expectation and application of the Gronwall lemma yield (|4.4I) . □ 

Corollary 4.3. The set {u^; e e (0,1)} is bounded in LP{n;C{[0,T];LP{T^))), 
for all p e [2, 00). 

Proof. To verify the claim, an uniform estimate of E(supg<j<2i ll''^'^(^)lliP(T«)) 
is needed. We repeat the approach from the preceding lemma, only for the 
stochastically forced term we apply the Burkholder-Davis-Gundy inequality. 
We have 

e( SUV \\u'(t)\\P„,^^A <E\\un\\P„,^^.+C{l+ I E\\u'(s)\\P„,^^As 



:(^sup^||u-(i)||i,(T«)) <IE||iio||^,(T«)+c(^l + E||u-(s)||P,(^„^ 

J2 f f \u'\P-\'gl{x,u')dxd(3k{s) 



+ p E sup 

\0<t<T 

and since we are dealing with finite-dimensional Wiener processes, the Burk- 
holder-Davis-Gundy, the assumption ()2.ip and the weighted Young inequality 
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yield 



sup 

0<t<T 



I I \u'r^u^9l{x,u-)<^xdPk{s) 

k=l -^T" 



k=l 

T . .2 

Lp(T«) 

T 



Therefore 

E(^sup^h-(OI|^,(T")) <c(^l+E|l^.o||^.(T«)+^'^Eh-(s)||^,(^„)ds 
and the corollary follows from ()4.4p . □ 
4.3. Compactness argument 

To show that there exists w : x x [0,r] — s- M, a kinetic solution to (|l.ip . 
one needs to verify the strong convergence of the approximate solutions . This 
can be done by combining tightness of their laws with the pathwise uniqueness, 
which was proved above. 

First, we need to prove a better spatial regularity of the approximate 
solutions. Towards this end, we introduce two seminorms describing the W^'^- 
regularity of a function u G L^(T^). Let A S (0, 1) and define 

Jtn JjN \x-y\^+^ 
Pg{u) = sup ^ / / \u{x) -u{y)\gr{x -y)dxdy, 

0<t<2Dn ''' JT" JT" 

where (gr) is a fixed regularizing kernel and by Dat we denote the diameter 
of [0, 1]-^. The fractional Sobolev space W^'^{T'^) is defined as a subspace of 
L^(T^) with finite norm 

l|w||M/Vi(TiV) = ll'uljii(Tiv) +p^{u). 

According to [6 , the following relations holds true between these seminorms. 
Let s E (0, A), there exists a constant C = Ca.^.at such that for all u G L^{T'^) 

[u) < Cp^ (u) , f<J--p^^(u). (4.6) 
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Theorem 4.4 (W^'^^^-regularity). Set a = min{^^^, i}, where a was introduced 
in (j2.3l) . Then there exists a constant Ct > such that for all t G [0, T] and all 
ee (0,1) 

Ep^(w^(t)) <Ct(1+Ep^K)). (4.7) 

In particular, for all s £ (0,cr), there exists a constant Ct,s,uo > such that for 
all te[0,T] 

E||u"(i)lllV-'.i(T") < CT,s,no(l+Eho|k-.i(T"))- (4-8) 

Proof. Proof of this statement is based on Proposition [XH We have 
E / / gr{x-y)nx,t,Of{y,t,Odidxdy 

<E [ [ Qr{x-y)M^-0r{x,t,0f{y,t,C)d^dCdxdy + 5 

<e/ / grix-y)M^-Ofo{x,Ofo{y,Od^dCdxdy + d + I + J + K 

<e/ / Qrix ~y)foix,Ofoiy,Od^dxdy + 2S + I + J + K. 

From the same estimates as the ones used in the proof of Theorem 13. 3[ we 
conclude 

e/ gr{x~y){u'{x,t)-u'{y,t))'^dxdy 

J(T")2 

<E / gr{x-y){uo{x)~uo{y))^ dxdy + 26 + Ct{6-^T + S^^T^ +6°') 

J(T")2 

By optimization in 6, i.e. setting S — , we obtain 

2S + Ct(S-W + S-^^ + S'') 
sup ^ < Ct, 

0<t<2Dn ^ 

where the maximal choice of the parameter a is min { , ^ } which corresponds 
to /? = max {7;^, 5}. Hence we deduce ()4.7p . Furthermore, by ()4.6p it holds 

Ep^(u^(t)) <Ct,.(1+Ep'^(7.o)) 
and from ()4.4p we obtain 

E||w^(t)|Ui(T«) < E||w-'(0||l2(t«) < Ct(i + (E||uo|li2(T«))^) 

and (14.81) follows. As a consequence of the previous estimate, the constant in 
(ji^ depends on the i^(ri; L2(T^))-norm of the initial condition. □ 

CoroUjiry 4.5. For all 7 £ (0, a) and q > 1 satisfying jq < a, there exists a 
constant C > such that for all e £ (0,1) 

Ell«'III.(o.T;W-.(T")) (4-9) 
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Proof. The claim is a consequence of the bounds (|4.4[) and (|4.8p . Indeed, fix 
s G (0, (t) and p G (1, oo). We will use the interpolation inequality 

<C\\-\\ (4.10) 

where 70,71 £ R-. Qo^Qi e (0,oo), 7 = (1 - 0)-fo + 9ji, 1 ^ l-± + 9 e {0, 1). 
In particular, fix s G (79,0') and set 70 = s, 71 = 0, qo = 1, gi = p. Then we 
obtain = , p = ^■''^^•''^ and 

(1-9)9 



<C(Eh-(i)l|v^=,i(T«)) (El|M-(i)||i,(T«)) <C- 

□ 

Also a better time regularity is needed. 

Lemma 4.6. Suppose that X G (0, 1/2), q G [2, 00). There exists a constant C > 
such that for all £ G (0, 1) 

''^ll^ llc'^([o,T];ff-2(ir")) — (4-11) 
Proof. Let q G [2,cxo). Recall that the set {u^; e G (0, 1)} is bounded in 

L«(f7;C(0,T; L?(T^))). 
Since all have the same polynomial growth we conclude, in particular, that 

{div(B^(M^))}, {div(A(a;)Vu")}, {eAu"} 
are bounded in ^^(fJ; C(0, T; iJ-^CjA^))) and consequently 



E 



9 



< c. 

Ci([0,T]:ff-2(T")) 



In order to deal with the stochastic integral, let us recall the definition 
of the Riemann-Liouville operator; let X be a Banach space, p G (l,oo], a G 
(1/p, 1] and / G iP(0,T;X), then we define 

{Raf){t) = ^^ J\t-sr-'f{s)ds, tG [0,r]. 

It is well known that Ra is a bounded linear operator from LP{0,T; X) to the 
space of Holder continuous functions C°'~^^p{[0,T]; X) (see e.g. [33 Theorem 
3.6]). Assume now that q G (2, 00), a G (1/(7,1/2). Then according to the 
stochastic Fubini theorem [5, Theorem 4.18] 

•P'{u'{s)) dW{s) - {RaZ){t), 



where 

r(i - a 



Z{s) ^ , / (s-r)-"r(u^(r))dW^(r). 
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Therefore using the Burkholder-Davis-Gundy and Young inequahty and the 
estimate (|2.ip 



9 



C°-i/5i([o,T];L2(T")) 
T , rt 1 



E 



< CTi^'-^^^Ej^ (l + h-(s)f^,(^„))ds 

< CTf (1-2") (l + h1IL(O;^,(0.T;L^(T«)))) < 

and the claim follows. □ 
CoroUjiry 4.7. For all d > there exist /3 > and C > suc/i t/iat /or all 

(0,1) 

IE||M"||c^([o,T];ff-''(T")) < C-. (4.12) 

Proof. The proof follows easily from interpolation between (14. lip and (14. 4p . □ 

CoroUjiry 4.8. Suppose that k G (0, 2(4*^0-) )■ There exists a constant C > smc/i 
i/iat /or aZ/ e e (0, 1) 

IEh1k»(0.T;L-(T«)) < C. (4.13) 



Proof. It follows from Lemma 14.61 that 



where A £ (0,1/2), q e [l,oo). Let 7 e (0, o-/2). If k = OX and = -26i+(l-6l)7 
then it follows by the interpolation and the Holder inequality 

IE||w''||if~(0,T;L2(T«)) < C'E(^||u^||^A(o_T;//-2(TiV))||w''||^2(o .r.^T(T«))) 

where the exponent r is chosen in order to satisfy (1 — 9)r = 2. The proof now 
follows from (|4J|) and (|4T4l) . □ 

Now, we have all in hand to show tightness of the collection {fi^; e £ (0, 1)} 
in X. Let us define the path space X = Xu x Xw, where 

Xu^L\0,T;L\T^))nC{[0,T];H-\T^)), Xw ^ C{[0,T];iX„) . (4.15) 

For all £ G (0, 1) we denote by the law of on X^ and by fiw the law of 
W on Xw- Their joint law on X is then denoted by fi'^ = ® iiw- 

Theorem 4.9. The set {/i"^; e G (0,1)} is tight and therefore relatively weakly 
compact in X. 



Degenerate Parabolic SPDEs 



23 



Proof. First, we employ an Aubin-Dubinskii type compact embedding theorem 
which, in our setting, reads (see (TS] for a general exposition; the proof of the 
following version can be found in [S]): 

L^{0,T;H'<{T^))nH''{0,T;L\T'^)) ^ L^{0,T;L^{T^)) ^ L\0,T; L\T^)). 
For i? > we define the set 

BiM = {ue L^{0, T; i7T(T^)) n H''(0, T; L^{T'^)); 

\\u\\l^{0,T-H-i{T")) + \\u\\h'^{O.T;L^{T^)) < R} 

which is thus compact in L^{Q,T; L\T'^)). Moreover, by gH) and 

f^u-{B?,R) < PM|w1lL2(0,T;//T(TiV)) > +PM|w''||h«(o,T;L2(T")) > ^ 



2 / \ C 



In order to prove tightness in C{[0,T]; H ^(T^)) we employ the compact em- 
bedding 

C^([0,T];i/-^(T^)) A C^{[0,T];H-\T^)) ^ C{[0,T]; H-\T^)), 
where 13 < (3, < -d < 1. Define 

B2,R ^{ue C^i[0,T];H-\T'')); \\u\\cfi^io,nH-^iT-)) < ^} 
then by (|iT^ 

f C 
fJ-u^{B2M) < ^1E||m''||c/3([o,t];H-''(t")) < 

Let ry > be given. Then, since Br — Bi r n B2.R is compact in and for 
some suitably chosen J? > it holds true 

Tj 

fJ.u-{BR) ^1^2' 

we obtain the tightness of {fiu^; e £ (0, 1)}. Since also the law fiw is tight as 
being a Radon measure on the Polish space Xw, there exists a compact set 
Crf C Xw such that ^w{C,j) > 1 — §• We conclude that Br x C,, is compact in 
X and fJ,^{BR x C^) > 1 — 77. Thus, {fj,^; e G (0, 1)} is tight in X and Prokhorov's 
theorem therefore implies that it is also relatively weakly compact. □ 

Passing to a weakly convergent subsequence /x" = (and denoting by ^ 
the limit law) we now apply the Skorokhod representation theorem to infer the 
following proposition. 

Proposition 4.10. There exists a probability space (fl,^,¥) with a sequence of 
X-valued random variables (u", VF"), n G N, and {u^W) such that 

(i) the laws of {ii"' ^W") and {u,W) under P coincide with /i" and 11, respec- 
tively, 

(ii) ({t",T4^") converges ¥-almost surely to {u,W) in the topology of X , 
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Note, that we can assume, without loss of generahty, that the cr-algebra 
^ is countably generated. This fact will be used later on for the application of 
the Banach-Alaoglu theorem. 

Remark 4.11. It should be noted that the energy estimates remain valid also 
for the candidate solution u. Indeed, let p € [2, oo) 

E sup < liminf E sup \\u"{t)\frp„N) 

0<t<T ^ ' 0<t<T ^ 

= liminfE sup \W{t)\\P„„. < C. 

n^oo o<t<T ^ 

Finally, let (=^f) be the P-augmented canonical filtration of the process 
(■u, W), that is 

.^t^<j{cT{gtU,QtW)u{N e.^: P(iV) = 0}), te [0,T], 
where gt denotes the operator of restriction to the interval [0,t\. 

4.4. Passage to the limit 

In this paragraph we provide the technical details of the identification of the 
limit process with a kinetic solution. The technique performed here will be used 
also in the proof of existence of a pathwise kinetic solution. 

Theorem 4.12. The triple (^t), P), T4^, u) is a martingale kinetic solu- 

tion to the problem (jl.ip . 

Let us define 

/" = 1u">Cj /" = 1«">^- / = lu>5i 
m" = < = (Vu")*yl(a;)(Vu")<5„.=e +e„|Vu"|^(5„.=j, 

m" = n1 + ^ (V{t")*yl(x)(Vt2")(5,v.=e + £„| Vw"|^(5s.=^. 

Let M.\y denote the space of bounded Borel measures over x [0,T] x R, i.e. 
the dual space of Ch^ the set of continuous bounded functions on x [0, T] x R. 

Lemma 4.13. It holds true (up to subsequences) 

(i) /"^/ in L°°{nxT^ x[0,T]xR)-weak*, 

(ii) there exists a kinetic measure rh such that 

fn^ — > fin in L^(n;Mb)-weak*. (4.16) 

Moreover, rh can be rewritten as hi+h2, where fii — (^u)* A{x){'S/ii)Su^^ 
and n2 is almost surely a nonnegative measure over x [0,r] x R. 

Proof. According to Proposition l4.10[ there exists a set S C 51 x x [0, T] of full 
measure and a subsequence still denoted by {u"; n e N} such that u^^uj, x, t) — > 
u{u!,x,t) for all {uj,x,t) e E. We infer that 

(4.17) 
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whenever 

(F^StCfN 0£[o,T]){('^,a;,i) e E; u{uj,x,t) = ^} = 0, 

where by C^n , C[o,t] we denoted the Lebesque measure on and [0,T], re- 
spectively. However, the set 



D 



e K; (f » Cj-N £[o,T]) {u^^)> o} 



is at most countable since we deal with finite measures. To obtain a contradic- 
tion, suppose that D is uncountable and denote 



N. 



Then D = UkenDk is a countable union so there exists /cq £ N such that Dfeg is 
uncountable. Hence 



I. °° 

fcn 

and the desired contradiction follows. We conclude that the convergence in 
(I4.17P holds true for a.e. (w, x, t, ^) and obtain (i) by the dominated convergence 
theorem. 

As the next step we shall show that the set {m"; n e N} is bounded in 
L^{(l; A4b)- Indeed, with regard to the computations used in proof of the energy 
inequality, we get from (|4.5I) 

Vu")*A{x){\/u")dxdt + en [ I |Vu"|^da;di < C||wo|li2(T«) 

-HCV / / u"g]i{x,u")dxdl3k{t) + C [ [ Gl{x,u")dxds. 
j^Jo Jt" Jo Jt" 

Taking square and expectation and finally by the Ito isometry, we deduce 

t|m"(T^ X [0,r] X R)|^ =E|m"(T^ x [0,T] x M)f 

/ / (VM")*A(a;)(Vw")da;dt + £„ / / \Vu"fdxdt 

Jo "'T" Jo "'T™ 



2 

< C. 



Thus, according to the Banach-Alaoglu theorem, (|4.16p is obtained (up to sub- 
sequence). However, it still remains to show that the weak* limit m is actually 
a kinetic measure. First, we conclude from (|4.4p 



t [ \^\P-^dfh"{x,t,0 < C 
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so for all fe > 

iin(|^|P~^fc)dm(x,i,C) 



rami 

T"x[0,T]xR 



= lim / min(|C|P-^fc)dm"(x,^,0 < C 

Jt"x[0,T]xR 

and consequently 

E [ dm"{x,t,0 < C. 

JT"x[0,T]xR 

Considering this fact and taking 0^ = -'-|C|>-R have by the dominated con- 
vergence theorem 

lim ]Em(T^ x [0,T] x B%) = E / lim (l3R{x,t,^)dm{x,t,^) = 0, 



where = G M; |^| > R}. Hence m vanishes for large ^. For predictability 
of 



JT«x[0,t]xE 

the same arguments as in can be used. 

Finally, by the same approach as above, we deduce that there exist kinetic 
measures 6i, 62 such that 

n" ^ 61, ^ 02 in L'^{i}; Mb)-weak* . 

Recall, that by <y{x) we denoted the square-root matrix of A{x). Then from 
()4.5p we obtain 

tf [ (VM")*A(a;)(Vu")dxdt = E /" [ \a{x){Vu")\'^dx dt < C 

Jo JT" Jo JT" 

hence application of the Banach-Alaoglu theorem yields that, up to subsequence, 
a{x){Vu") converges weakly in L^{(l x x [0,T]). On the other hand, from 
the strong convergence given by ProDOsition l4.10l we conclude using integration 
by parts, for all V e C^iT^ x [0,r]), 

a{x){\/u")i:{x,t)dxdt — >[ [ a{x){Vu)tp{x,t) dx dt, P-a.s.. 



10 JT" Jo JT" 

Therefore 

(7(x)(V7-i") cr{x){\/u), in L^{T^ x [0,r]), P-a.s.. 

Since any norm is weakly sequentially lower semicontinuous, it follows for all 
<fi e Cb(T^ X [0,T] X R) and fixed ^ G M, P-a.s., 

/ / \a{x){'^u)\'^^^{x,t,^)dxdt <\imM [ [ \a{x)i\7u")\^ip^{x,t,^)dxdt 

Jo Jt" n-s-oo Jq Jjn 
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and by the Fatou lemma 

Iff \(j{x){\/u)\^ip'^{x,t,£_)dSu=^dxdt 
Jo Jt" Js. 

<liminf/ / [ \a{x){Vu^')\'^(p^{x,t,£_)dSu^^^dxdt, P-a.s.. 
Jo Jt" Jm 

In other words, this gives ni = {\/u)* A{x){\/u)6u=^ < oi P-a.s. hence h2 = 
02 + (oi — is P-a.s. a nonnegative measure and the proof is complete. □ 

Let us define for all t e [0, T] and some fixed if G C^{T^ x M) 

M"(t) = (rw, ^) - (/o, - /* (r(s), 6"(e)- v^) ds 

Jo 

ft , ^ , /-t 

-/ (r(s), 5] a,,(A,,(x)a.,^))ds-£„ / {r{s),Aip)ds 

Jo \ ,^^-^1 / Jo 

{Sw^^^Gl, d^ip) ds + {m'\d^ip)i[0, t)), n e N, 

Jo 

- / (/"(^)' E 5.,-(Ay(x)a,.(p))ds-e„ / (r(s),A^)d5 
Jo \ ,,^^.^1 / Jo 

(5«"=«G2 , d^ip) ds + {m'\d^ip)i[0, t)), nen, 

M{t) = (fit), cp) - (/o, cp) - f {f{s)MO- V^) ds 

Jo 

Jo ^j^^ Jo 

+ (m,a5^)([0,i)). 

The proof of Theorem 14.121 is an immediate consequence of the following two 
propositions. 

Proposition 4.14. The process W is a (.^t)- cylindrical Wiener process, i.e. there 
exists a collection of mutually independent real-valued (^t) -Wiener processes 
{i3k}k>i such that W = J2k>i Pkek- 

Proof. Hereafter, times s,t g [OiT^li s < t, and a continuous function 

7:C([0,s];i/-i(T^)) xC7([0,s];ilo) ^[0,1] 

will be fixed but otherwise arbitrary and by we denote the operator of re- 
striction to the interval [0, s]. 
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Obviously, 1^ is a ilo-valued cylindrical Wiener process and is (^t)-adap- 
ted. According to the Levy martingale characterization theorem, it remains to 
show that it is also a (^t)-martingale. It holds true 

Ej{gsu'\ QsW^) [W^{t) - VK"(s)] - E^{g,u'\ g,W) [W{t) - W{s)] = 

since is a martingale and the laws of (w", W") and (u", W) coincide. Next, 
the uniform estimate 

sup]E||W^"(i)||i - supE||W^(i)||i < oo 

nGN neN 

and the Vitali convergence theorem yields 

Ej{gsu, g,W) [W{t) - W{s)] = 
which finishes the proof. 

Proposition 4.15. The processes 



□ 



/ {Su^^gk,ipy dr, M/3fc - / {St,^^ gk,(p) dr 



are (.^t)-'rncLrtingales. 

Proof. Here, we use the same approach and notation as the one used in the 
previous lemma. Let us denote by fc > 1 the real- valued Wiener processes 
corresponding to W^", that is = X]fc>i f^k^k- For all n G N, the process 

M"- / (5„n=^<?"K)dW^,^) / (<5„..^g,",^)d/3fc(r) 

Jo f^^j^Jo 

is a square integrable (^i)-martingale by (|2.ip and by the fact that the set 
{m"; n e N} is bounded in L^{n; L^{0,T; L^{T^))). Therefore 



(M")2 f g^, dr, M"/3fe - / (<5„.=e 9k, v) 



dr 



are (^i)-martingales and this implies together with the equality of laws 

t^{g,u'\ g,W^) [M"(0 - M"(s)] = E-i{g,u'\ g,W) [Nf^t) ~ Af"(s)] = 0, 

(4.18) 



{M-fit) - {£rf{s) -Y^f (5,.=^ 9l, dr 
(M")2(t) - (M«)2(s) - ^ r (^„„^^ g^, dr 



fe>i 



= 0, 
(4.19) 
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= 0. 
(4.20) 

Moreover, the expectations in (|4.18p - (|4.20p converge by the Vitali convergence 
theorem. Indeed, all terms are uniformly integrable by (|2.ip and ()4.4p and con- 
verge P-a.s. (after extracting a subsequence) due to Lemma I4.13| Proposition 
14.101 and the construction of , . Hence 



t-f{gsU,gsW)[M{t)-M{s)\ = 0, 



E-^{gsU,gsW) 



E-f[gsU,gsW) 



k>l''^ 



= 0, 



M(i)/3fe(t) - M{s)pk{s) - / {6u=igk,^)dr 



= 0, 



which gives the (^4)-martingale property. 



□ 



Proof of Theorem 14.121 Once the Propositions 14.141 14.15l are established, it fol- 
lows that the quadratic variation 



M 



and so for every (p e Cf=(T^ xR), t£ [0,T], P-a.s 



(/>),^)-(/o,¥')- / (/>),6(e)-V^)ds- / (/>), ^a.^.(A,,(x)9.,^))ds 



= / {Su^^'P{u)dW,^) + ^ I {6i,=^G^d^^)ds-{rh,d^^){[0,t)) 



and the statement follows. 



□ 



4.5. Pathwise solutions 

In order to finish the proof, we make use of the Gyongy-Krylov characterization 
of convergence in probability introduced in |12j . It is useful in situations when 
the pathwise uniqueness and the existence of at least one martingale solution 
imply the existence of a unique pathwise solution. 

Proposition 4.16. Let X be a Polish space equipped with the Borel a-algebra. A 
sequence of X -valued random variables {Yn; n G N} converges in probability if 
and only if for every subsequence of joint laws, {fJ-nk,mkj k G N}, there exists a 
further subsequence which converges weakly to a probability measure fx such that 



fi{{x,y) eX xX;x = y) =1. 
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We consider the collection of joint laws of {u",u"^), denoted by /J.J^'™. For 
this purpose we define the extended phase space (cf. (|4.15l) ') 

X'^ = Xu X Xu X Xw, X^ = Xu X Xu^ 

As above, denote by /i" the law of u" and by jiw the law of W . Set further 

Similarly to Proposition 14.91 the following fact holds true. The proof is 
nearly identical and so will be left to the reader. 

Proposition 4.17. The collection {;/"'"'; n, m G N} is tight on X K 

Let us take any subsequence {^"'='™fe; fc £ N}. By the Prokhorov theorem, 
it is relatively weakly compact hence it contains a weakly convergent subse- 
quence. Without loss of generality we may assume that the original sequence 
jj^Mfc, "!■):• g f^} itself converges weakly to a measure v. According to the Sko- 
rokhod representation theorem, we infer the existence of a probability space 
($1,^,1') with a sequence of random variables (u"*" , u™'' , T^'^), fc e N, conver- 
ging almost surely in X^ to a random variable (u, u, W) and 

P((^i"^^i"^T^'=) e •) = j/"'''™^-), f{{u,u,W) e •) = 
Observe that in particular, converges weakly to a measure /i„ defined 

by 

^iu{■) = f{{u,u) e •)■ 

As the next step, we should recall the technique established in the previ- 
ous section. Analogously, it can be applied to both , 14^*^), {u^W) and 
W^), {H, W) in order to show that (u, W) and (u, W) are martingale ki- 
netic solutions of (jl.ip defined on the same stochastic basis (f2, (^t), P), 
where 

^ cT{c7{gtU, gtii, QtW) U {N e ¥{N)^0}), te [0,T]. 

Since u(0) = ii{0) P-a.s., we conclude from Theorem [33] that u = iiin Xu P-a.s. 
hence 

lJ-u{{x,y) e Xu X Xu] X = y) = P(u = u in Xu) = 1. 

Now, we have all in hand to apply Proposition 14.161 It implies that the origi- 
nal sequence u" defined on the initial probability space {il,J^,P) converges in 
probability in the topology of Xu to a random variable u. Without loss of gene- 
rality, we assume that m" converges to u almost surely in Xu and again by the 
method from section we finally deduce that u is a pathwise kinetic solution 
to (|l.ip . Actually, identification of the limit is more straightforward here since 
in this case all the work is done for the initial setting and only one fixed driving 
Wiener process W is considered. 
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5. Existence - general initial data 

In this final section we provide an existence proof in the general case of uq G 
Z/P(f2; LP(T^)), for all p £ [l,oo). It is a straightforward consequence of the 
previous section. We approximate the initial condition by a sequence {uq} C 
LP{n;C°^{T^)), p G [l,oo), such that ^ uq in (ft; (T^ )) . That is, the 
initial condition Uq can be defined as a pathwise mollification of uq so that it 
holds true 

IWoW LP {n; LP iT«)) < ll"o||LP(n;LP(T«)), £ G (0, 1), P G [1, oo). (5.1) 

According to the previous section, for each e G (0,1), there exists a kinetic 
solution to (jl.ip with initial condition Uq. By application of the comparison 
principle 



E\\u'' -u'-Ui^o,T;L^iTN)) <TE\\u'^' -uI'Wlut-), £i,£2 G (0,1). 

Therefore, {u^; e G (0, 1)} is a Cauchy sequence in L^{n; ^^(O, T; ^^(T^))) and 
there exists u G L^{Vt;L^{{),T]L^{T'^))) such that 

— >u in Li(17;Li(0,r;i^(T^))). 

By ()5.ip . we still have the uniform energy estimates 

E sup \\u%t)\Wj,,^.<CT,uo, pG[2,cx}), /t-o\ 

0<t<T ^ ' ^ ' 

as well as (using the usual notation) 

E|m'^(T^ X [0,r] X R)|^ < Ct,uo- (5.3) 

Thus, using this observations as in Lemma I4.13| one finds that there exists a 
subsequence {-u"; n G N} such that 

(i) /" ^ / in i°°(l] X X [0,T] x R)-weak*, 

(ii) there exists a kinetic measure m such that 

m" ^ m in i^(f2; A4(,)-weak* 
and m = ni+n2, where ni = {\7u)* A{x){'S/u)du=(^ and n2 is almost surely 



a nonnegative measure over T x [0, T 



X 



With these facts in hand, we are ready to pass to the limit in (|2.7p and conclude 
that u satisfies the kinetic formulation in the sense of distributions. Note, that 
(|5.2p remains valid also for u so (|2.6p follows and the proof of Theorem 12.91 is 
complete. 
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